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AN  APPROXIMATE  METHOD  OF  CALCULATING  DOWN\v’ASH  BEHIND  A STRAIGHT  WING  WITH 
UNSTEADY  APERIODIC  MOTION  AT  SUBSONIC  FLIGHT  SPEEDS 

0.  M.  Danchenko 

Information  about  downwash  angles  behind  a wing  of  the  specified  type  and 
the  character  of  their  relationship  to  the  angle  of  attack  of  the  wing  during 
steady  and  unsteady  motion  is  vital  for  solving  problems  of  the  efficient  con- 
figuration of  the  tail  unit. 

The  calculation  of  downwash  at  the  aircraft  tail  unit  is  a comparatively 
simple  problem  under  stationary  conditions  Dl.  Unsteady  downwashes  created 

A. 

by  a wing  are  more  difficult  to  determine. 

The  investigation  of  downwash  created  by  a wing  during  harmonic  vibration 
was  presented  in  a quite  general  way  in  a work  [23,.  The  question  of  investi- 
gating unsteady  downwashes  during  aperiodic  motion  is  of  equal  practical  interest 
because  one  must  study  the  motion  of  an  object  with  aperiodic  relationships  of 
the  kinematic  parameters  and  time  when  the  aircraft  vertically  zoom-climbs 
(breaks  away),  changes  from  one  angle  of  attack  to  another,  etc. 

This  article  presents  a numerical  method  of  calculating  downwash  angles  at 
the  aircraft  tail  unit  for  unsteady  aperiodic  motion  at  subsonic  speeds.  Only 
straight,  untwisted  wings  are  examined,  but  it  seems  possible  to  devise  a similar 
method  for  wings  of  any  shape  with  aerodynamic  or  geometric  twist. 

A coupled  system  of  rectangular,  rectilinear  coordinates  xyz  is  introduced 


in  which  the  x axis  runs  backward  along  the  medial  aerodynamic  chord  of  the 

wing,  the  z axis  is  to  the  right,  along  the  wing  span  on  a line  perpendicular  to 

the  plane  of  symmetry  of  the  wing,  and  the  y axis  is  above,  perpendicular  to  the 

xz  plane.  The  origin  of  the  coordinates  is  at  ^ CAX  (Figure  1,  a).  Furthermore, 

4 

we  shall  use  a relative  point-grid  reference,  having  taken  half  the  wing  span 
as  the  typical  dimension: 


We  replace  the  straight  wing  with  a system  of  bound  vortices  with  axes 

parallel  to  the  z axis.  These  vortices  traverse  the  pressure  center,  whose 

location  we  assume  to  be  i CAX;  their  combined  circulation  is  equal  to  true 

4 

circulation  in  each  wing  section.  A sheet  of  free  vortices  whose  axes  are 


parallel  to  the  velocity  of  the  free-stream  flow  if  viewed  from  above  streams 


I 

! 


n L' 


away  from  the  wing.  This  sheet  is  unstable  and  turns  a certain  distance  behind 
the  wing  into  two  vortex  cores  (Figure  1,  a)  as  the  result  of  the  mutual  effect 
of  the  vortices  that  form  it.  One  can  assume  that  the  process  of  turning  is 
complete  in  the  vicinity  of  the  tail  unit  and  that  the  vortex  cores  begin 
immediately  adjacent  to  the  wing  [l,  Sj  . Such  an  arbitrary  vortex  model  sig- 
nificantly simplifies  solving  the  problem  and  is  only  valid  at  an  adequate  dis- 
tance from  the  wing,  which  exists  in  the  case  of  tail  unit  configuration  on  the 
aircraft. 

If  the  vortex  circulation  does  not  change  in  time,  then  the  given  vortex 
model  corresponds  to  steady  motion.  In  this  case  the  distance  between  the  free 
vortices  depends  on  the  distribution  law  of  circulation  along  the  wing  span  and 
is  calculated  on  the  basis  of  the  N.  Ye.  Zhukovskiy  theorem,  according  to  which 
lifts  determined  near  the  wing  and  far  behind  the  wing  where  the  effect  of  the 


bound  vortices  can  be  ignored  are  identical: 


(1.1) 


whence 
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where  1 is  wing  span; 

1^  = vl  is  the  distance  between  free  vortices; 

r = f(z)  is  the  law  of  distribution  of  velocity  circulation  along  the  wing  span; 
Tq  is  velocity  circulation  in  the  plane  of  symmetry  of  the  wing. 

The  law  of  circulation  distribution  along  the  span  of  a straight  wing  is 
assumed  to  be  the  average  between  trapezoids  (with  respect  to  the  base  lengths, 
equal  to  k - a value  which  is  the  reciprocal  of  the  (wing)  taper  ratio  and  an 


ellipse : 


r = o.5rol]/i-c»+  i-(i-y^)C), 


(1.2) 


k = 1 - for  a parallel  wing; 
k = 0 - for  a delta  wing. 

According  to  the  calculations  cn  , the  accepted  law  of  circulation  distri- 
bution coincides  with  the  true  law  with  sufficient  accuracy.  One  can  express  cir- 
culation in  the  plane  of  symmetry  of  a wing  through  the  lift  coefficient  and 


With  the  selected  law  of  wing  span  circulation  distribution,  the  distance 
between  free  vortices  depends  solely  on  the  geometrical  characteristics  of  the 
wing.  Actually,  on  the  basis  of  (1.1),  (1.3),  and  (l.iJ) 

1^  = 0.51  (1.285  + 0.5k).  (1. 

Circulation  of  the  bound  vortices  is  constant  in  accordance  with  the 
constant  values  of  the  angle  of  attack  and  angular  velocity  during  steady 
motion.  It  changes  in  time  during  unsteady  motion  and  this  can  be  inter- 
pretted  as  the  result  of  the  formation  of  an  additional  rectilinear  closed 
vortex  line  with  length-wise  uniform  circulation  AI  (Figure  1,  b).  The  trans- 
verse vortex  with  circulation  - AT  - is  carried  downstream  and  the  distance 
from  it  to  the  bound  vortices  constantly  increases,  i.e.,  a free  vortex  of  a 
new  type  - a transient  free  vortex  - appears.  Its  downwash  with  a velocity  V 
in  the  plane  formed  by  the  free  steady  vortices  changes  the  circulation  of 
the  latter  length-wise  and  in  time. 

In  order  to  solve  the  problem  in  the  case  of  an  unsteady  aperiodic  motion, 
we  shall  present  the  continuous  process  of  change  of  the  aerodynamic  character- 
istics of  the  wing  in  time  in  the  form  of  a set  of  discrete  changes.  The  kine- 
matic parameter,  and  consequently,  the  aerodynamic  characteristics  change  jump- 
wise  at  certain  moments  of  time  and  the  intervals  between  these  moments  are 
constant  (Figure  2). 

We  choose  the  calculated  moments  of  time  (t^,  t^^,  ...,  t^)  such  that  they 
directly  precede  the  moments  in  which  the  kinematic  parameters  change  jump-wise, 
and  thus,  the  circulation  of  the  bound  vortices. 

At  the  beginning  of  unsteady  motion  (t^),  the  circulation  of  vortices  is 
determined  by  steady  initial  values  of  the  kinematic  parameters  and  the  vortex 
system  has  the  form  shown  in  Figure  1,  a.  We  shall  call  it  a quasi-steady 
vortex  system. 

Inasmuch  as  t^  is  the  beginning  of  unsteady  motion,  then  there  are  no 
transient  vortices  at  this  calculation  moment  and 

c*,®' = -t- (1.6) 

where  and  are  the  initial  values  of  the  kinematic  parameters  which 

3 z 

correspond  to  the  beginning  of  unsteady  motion;  the  parenthetical  0 as  a super- 
script signifies  the  calculated  moment  of  time  when  the  value  of  the  corres- 
ponding parameter  is  determined: 
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Figure  2. 


Circulation  of  the  bound  vortices  will  change  suddenly  by  a moment  of  time 
t^,  which  will  be  accompanied  by  convergence  of  the  transverse  free  vortex  that 
will  recede  a distance  VAt  from  the  bound  vortices  by  this  time. 

The  lift  factor  at  a calculated  moment  of  time  t^  is 


where  and  values  of  the  kinematic  parameters  after  a single 

sudden  change  that  corresponds  to  a moment  of  time  t^. 

On  the  one  hand,  the  velocity  induced  at  any  point  can  be  viewed  as  the 
result  of  the  effect  of  the  quasi-steady  vortex  system  that  corresponds  to  a 
calculated  moment  of  time  t^  and  the  additional  closed  cortex  system  shown  by 
the  dash  line  in  Figure  1,  b.  But  on  the  other  hand,  which  is  more  convenient, 
one  can  consider  that  a new  quasi-steady  vortex  system  has  formed  at  a calculated 
moment  t^  which  has  the  previous  law  of  distribution  of  circulation  of  the  bound 
vortices  along  wing  span  given  by  expression  (1.2)  with  a new  value  of  the  lift 
factor. 

Two  free  vortex  cores  with  uniform  circulation  + AT^  run  to 

infinity  from  the  bound  vortices,  but  then  the  transverse  free  circulation 


/lO 
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vortex  - AP^  - also  has  longitudinal  free  vortex  lines  that  also  run  to  infinity, 
i.e.,  we  obtain  a new  quasi-steady  vortex  system  and  one  horseshoe  vortex  with  a 
length  (span)  1^  (Figure  1,  c). 

The  magnitude  of  the  circulation  change  AT^  can  be  calculated  in  the  follow- 


ing way: 


if. ri"  _ ri"  » ^ _ Ar„ 


(1.8) 


The  above  is  also  valid  when  there  is  not  one  but  a whole  series  of  transient 
transverse  free  vortices  behind  the  wing.  At  each  succeeding  calculated  moment 
of  time  t^,  the  horseshoe  vortices  examined  at  the  previous  moment  of  time  t^  -1 
move  downstream  a distance  VAt  without  a circulation  change,  and  moreover,  a 
new  horseshoe  vortex  appears  with  a circulation  AF^.  The  horseshoe  vortex 
corresponds  to  the  r-th  change  of  the  quasi-steady  vortex  system  and  will  cover 
a distance  VAt  from  the  wing. 

n changes  of  the  quasi-steady  vortex  system  occur  and  n systems  of  horseshoe 
vortices  forms  by  a moment  of  time  t . At  a moment  of  time  t , the  kinematic 
parameters  acquire  values  of  and  tii'  \ and  with  this 

3 Z 


= 4 . 4 r', . 

o t 2/j  ’ 

Aa,  = — a','’-'!;  AuT,„  = •Hi'”  — 


Ar„  = Ac, 


2/j  • 


The  given  vortex  model  makes  it  possible  to  calculate  the  downwash  angle 
behind  the  wing  during  unsteady  aperiodic  wing  motion, 

2.  We  shall  solve  a problem  in  the  linear  presentation,  i.e.,  we  shall 
consider  the  sheet  of  free  vortices  behind  the  wing  to  be  plane  and  the  free 
vortex  cores,  which  represent  the  results  of  twisting  of  the  vortex  sheet,  to 
be  rectilinear  and  lying  the  the  xOz  plane.  The  linear  presentation  signifi- 
cantly simplifies  the  problem  and  does  not  result  in  large  errors  in  the  region 
of  low  angles  of  attack  of  the  wing  hJ  . This  region  is  of  the  greatest  prac- 
tical interest,  for  the  low  angles  of  attack  of  the  wing  are  the  critical  ones 
for  the  stabilizer-elevator  unit  (in  the  sense  of  stall  probability). 


- 6 - 


We  shall  consider  that  the  wing  moves  translationally  at  a uniform  velocity 
Vq,  that  flow  braking  behind  the  wing  is  aosent  and  that  the  downwashes  created 
by  the  wing  are  constant  along  the  tail  span.  We  determine  the  projections  of 
velocity  in  the  direction  perpendicular  to  the  free-stream  flow  (direction  mm). 

These  projections  are  induced  by  elements  of  the  vortex  model  at  calculation 
point  A (C,  n,  0),  which  is  located  in  the  plane  of  symmetry  of  the  wing  (Figure  3). 

On  the  basis  of  a work  1 , the  velocity  induced  by  the  system  of  bound 

vortices  distributed  along  the  wing  span  according  to  a law  (1.2),  is: 

[l*  + l+i>+,.-KF+T>  . 


Figure  3. 

Substituting  (i.U)  in  this  expression,  we  obtain 

‘^^0  > 1 r , 


■ 2r.\  ('1.285  + 0,5*) 


1 

nwr 


+ T,* 


Where 


- ( 1 - A:)  KTTF+v  - . 

a 

y is  a function  of  the  total  elliptical  integrals  of  the  first  and  second 


orders . 


The  projection  of  velocity  to  direction  mm  is 


“'io-=U',«-cos(a  — Po). 


(2.2) 
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where 


6.  = arctg  _ . 

C 

Angle  £ is  nearly  zero  for  aircraft  whose  stabilizer-elevator  units  are  on 


the  fuselage.  In  this  case 


uy,o  = aKioCosa. 


In  the  region  of  small  angles  of  attack  of  the  wing 

CO  = OJ  . 

10  10 

Velocity  induced  by  two  vortex  cores  with  circulation  1^  is 


(2.3) 


(2.M 


where 


1 


TV  ''  * 

Bearing  in  mind  that  ^0=  — . we  express  oi  through  coefficient  Cy. 


2nx  v«  4-  T.*  ( ' + y^r:^  ■ 


(2.5) 


T'  •■’Ction  of  velocity  to  direction  mm  is 

“20  = “20 

and  in  t-he  region  of  small  angles  of  attack  is 


(2.6) 


“20  ■ “20' 


(2.7) 


Velocity  induced  by  a free  transverse  vortex  with  circulation  - Air,  is 

-ai,  1 


where 


Taking  expression  (1.8)  into  account,  we  write 

_ ' 1 

2«.>  y V/v>  + e»  + T*‘ 

The  projection  of  velocity  to  direction  mm  is 
tti„  = wi,  • cos  (a  — cosa  • cosp,. 


where 


p,  = arctg  . 
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(2.8) 


(2.9) 


L, 


(2.10) 


In  the  region  of  small  angles  of  attack 

= tt'u  cos  p,  = W\,  - . . 

K + t* 

Velocity  induced  by  two  longitudinal  vortices  with  circulation  - Afr  is 


nt  + T*  ^ j • 


Substituting  expression  (1.8),  we  have 

Ac..  V. 

tVir 


, = ( 1 4- ■ - 

2«A  + }/>«  + {»  4 


The  projection  of  velocity  to  direction  irnn 


The  region  of  small  angles  of  attack 

102,  = W2r  COS  a. 

Rake  angle  at  a calculation  point  A is 

Ws,  = Wg,. 

because  of  the  effect  of  the  vortex  systems  exa.mined  above,  where 
(1.285  + 0.5*)  |'c“+>,*  +^*+V  — 


J + v*+  ^ 4.  j 


y.r  = 


ir 

+ > 

f.  , ] 

''  (-  + + r,s 

' >*  + V 

L ' Kv*  + e2  4.,.J 

In  the  region  of  small  angles  of  attack 

^KP  = i^(^pXo  — Ac^.Xr)- 


(2.11) 

(2.12) 

(2.13) 

(2. ID) 


(2.15) 

(2.16) 


(2.17) 


Value  Xq  is  solely  a function  of  the  coordinates  of  the  calculated  point 

and  does  not  depend  on  time,  and  x^,  is  a function  of  the  coordinates  of  the 

calculated  point  and  value  C^.  With  assigned  wing  velocity  and  a chosen  span 

of  time  At,  value  C changes  according  to  a linear  law  that  multiply  corresponds 
r 

to  calculated  moment  r (r  = 0,  1,  2,  ...,  n). 

Taking  the  above  into  account,  it  is  easy  to  devise  a method  of  calculating 
the  downwash  angle  behind  a wing  in  the  region  of  the  stabilizer-elevator  unit 
during  unsteady  aperiodic  motion.  Actually,  by  successively  examining  the 
corresponding  vortex  models,  we  calculate  the  projection  of  velocity  induced 
by  all  bound  and  free  vortices  at  each  calculated  moment  of  time  to  direction  mm. 


and  then  also  calculate  the  corresponding  downwash  angle. 

The  vortex  nodel  appears  as  shown  in  Figure  1,  a at  a calculated  moment 
of  time  t^  (r  = 0)  at  the  beginning  of  unsteady  motion,  i.e.,  it  consists  solely 
of  a quasi-steady  vortex  system.  Designating  the  calculated  moment  of  time  with 
the  index  in  parentheses  and  taking  into  account  that  = 0 when  r = 0,  one 

can  write  the  following  on  the  basis  of  (2.14) 

,,0)  _ + ‘-’20  _ y 

cr  ~ 

where  is  calculated  for  values  of  the  kinematic  parameters  and 

_(i)  y a 

, and  calculated  according  to  formula  (2.16)  when  r = 1; 

Ac„.  - c'  — 


The  vortex  model  consists  of  a quasi-steady  vortex  system  and  horseshoe 

vortices  whose  number  is  equal  to  p at  a calculated  moment  of  time  t . 

P 

Downwash  angle 

t<p>  = . -+  + S'/;  + . . . + 

cr  V'o  ~ ~ ~ = 

~ 27./.  IV  Xo  — Ac^,Xp  — ACp,Xp_,  — . . . _ _ . . . _ Arp^Xi]  = 

= V J 

2r.A  ‘-p  '■»  L I • 

1=1 


(2. 17) 


In  the  region  of  small  angles  of  attack  of  the  wing 


e'p)  = _L 
cr 


kWxo-  V V Ac^y,'. 

!=\  i^.p 


(2.18) 


where 


4'”  = /IV’.  V’]; 

ACp,  = cy)_c(0); 

ACpp  = C(r)^f(^l); 

The  values  x-j^  > X2>  •••>  X calculated  according  to  formula  (2.16)  with 

the  corresponding  values  r = 1,  2,  ...,  p. 

3.  At  high  subsonic  flight  velocities,  one  can  calculate  downwash  angles 
behind  the  wing  using  the  Prandtl-Glauert  transform. 

We  transform  the  coordinate  system  such  that 


xm  = -7—  • ■ : hm  — y<  2m  — 2. 

K 1 — Mi 


(3.1) 
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I 


where  x , y.,  and  z are  the  coordinates  of  a point  in  transformed  space  which 
M ” M 

corresponds  to  the  calculated  one. 

The  geometrical  parameters  of  the  wing  in  the  transformed  space  are: 

^ ]/l  — M^:  ca(  =c  ; Im  = 1\  bs\  = - ■■  ; = — =^=r.  (3  2 

We  shall  take  the  following  kinematic  parameters  for  the  transformed  wing: 


1/  ..Ur 

OAf  — »'  0*  — a; 


(here  “’^a)  ~ ) > then  the  inductive  velocities  in  the  corresponding 

points  of  the  original  and  transformed  spaces  are  linked  by  relationships  Csl 


I I — A\  J 


(3.3) 


and  the  following  relationships  exist 


^VA.  = = c?  y 1 - Mi;  -f  Cv 


(3.'4) 


We  calculate  the  parameters  of  the  transformed  wing  and  the  coordinates  of 
the  corresponding  point  according  to  the  assigned  geometrical  parameters  of  the 
wing  and  the  coordinates  of  the  calculated  point  according  to  formulas  (3.2)  and 
(3.1).  Then  we  determine  the  inductive  velocities  in  transformed  space  on  the 
basis  of  expressions  (2.1),  (2.5),  (2.8),  and  (2.11).  Knowing  the  inductive 
velocities  according  to  formulas  (3-3)  and  (2.17),  we  find  the  downwash  angles 


behind  the  actual  wing  at  the  given  Mach  number. 


E<P>  = 

cr 


. - Si  [ ‘Zi..  - 1 i J - 


(3.5) 


4.  The  angle  of  attack  of  the  stabilizer-elevator  unit  in  steady  curvilinear 

flight  is  determined  by  the  following  relationship: 

a =a  + (()-c  -E„  + Aa  , (4.1) 

s . c . cr  0 s . c. 


where  a is  the  angle  of  attack  of  the  wing; 

is  the  rigging  angle  of  incidence  of  the  stat ilizer-elevator  unit; 

G is  the  downwash  angle  caused  by  the  wing; 
cr 


e„  is  the  downwash  angle  when  C =0. 

0 y 


*T  --V  .. 
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Figure  4. 


Aa,.o=^'o$- 


- the  additional  change  of  the  angle  of  attack  of  the  stabilizer-elevator  unit 

caused  by  turning  of  an  aircraft  with  angular  velocity  w . 

z 

One  can  calculate  the  angle  of  attack  of  the  stabilizer-elevator  unit  by 
the  numerical  method  during  unsteady  aperiodic  motion  in  the  vertical  plane 
using  formulas  (4.1)  and  (4.2),  having  replaced  the  continuous  processes  with 
discrete  ones.  Successively  examining  the  vortex  models  that  correspond  to 
calculated  moments  of  time  t^,  t^,  ...,  t^,  we  determine  the  downwash  angles 
^(0)^  ^(1)  ^(p)  additional  changes  of  the  angle  of  attack  of  the 

stabilizer-elevator  unit  Aa^  Aa^  Aa_  according  to  formula  (2.17) . 


The  angle  of  attack  of  the  stabilizer-elevator  unit  at  a calculated  rjonent 


of  time  t is 
o 


where 


au»  = oio»  + <?  — — £,  + A»v»)  , 
s.e  cr  s.e 


= a<»»  — a,; 


s.e  K i '■ 

Figures  4-6  show  oscillograr.s  with  recordings  if  v 
of  aircraft  motion  when  the  aircraft  is  maneuvering  ir, 
it"  the  elevator).  The  operating  conditions  of  tne  e'^f.r, 
(B  = 0);  flaps  up. 
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Figure  5. 


r 


The  coordinates  of  the  calculated  point  are:  K = 0.697;  n = 0.093;  C = 0. 

The  wing  parameters  are:  X = 12;  k = 0.3814;  = 1°. 

The  results  of  calculations  made  according  to  the  suggested  method  are  in 
good  agreement  with  the  data  of  a flight  experiment  carried  out  in  broad  ranges 
of  angle  of  attack  and  Mach  number. 


Figure  6. 
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